Quantum Theory of Electrodynamics in Linear Media Subject to Boundary 

Conditions 
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We show that the material dependencies of macroscopically quantized fields in linear media are 
not consistent with the classical electromagnetic boundary conditions. We then phenomenologically 
construct macroscopic quantized fields that satisfy quantum-classical correspondence with the re- 
sult indicating that the canonical momentum in a linear medium is modified as a consequence of 
the reduced speed of light. We re-derive D'Alembert's principle and Lagrange's equations for an 
arbitrarily large region of space in which light signals travel slower than in the vacuum and show 
that the resulting modifications to Lagrangian dynamics, including the canonical momenta, repair 
the violation of the correspondence principle for macroscopically quantized fields. 



I. INTRODUCTION 

Quantum electrodynamics describes the quantum me- 
chanical properties of electromagnetic fields in the vac- 
uum. The analogous theory for the quantum optical 
properties of electric and magnetic fields in a dielec- 
tric or similar linear medium can be obtained by quan- 
tizing the macroscopic fields in the continuous medium 
P, H H, 0, H, as well as by applying quantum 
electrodynamics to a microscopic model of the medium 
0, S, B El- Adopting the former approach in 1940, 
Ginzburg [l| quantized the macroscopic electromagnetic 
field in a simple transparent linear dielectric. The ide- 
alized dielectric studied by Ginzburg is a model for a 
more complex material that exhibits dispersion, absorp- 
tion, nonlinearity, and other extended effects to vary- 
ing degrees. The difficulties of treating dispersion within 
the macroscopic quantization procedure have been ad- 
dressed by Jauch and Watson Drummond [s^l, and 
Milonni [4j|. Huttner and Barnett ^] treated absorptive 
dielectrics using a macroscopic quantization procedure 
based on damped Hopfield polaritons. Quantization of 
the macroscopic field in nonlinear dielectrics [Tl| , magne- 
todielectric matter |12|], and negative-index materials 13 1 
has also been investigated. In contrast to reports of ex- 
tended aspects of quantized fields in matter, the current 
work challenges the validity of the macroscopic quanti- 
zation procedure itself. Specifically, we use the classical 
electromagnetic boundary conditions to show that the 
quantized field inside a linear medium, as it is currently 
derived P, d, ^, ^, ^, 6.] , violates quantum-classical cor- 
respondence. 

Materials are of finite extent and we are usually in- 
terested in how the dynamics of a field change with the 
properties of the medium. In the limit of large num- 
bers, the quantum and classical relations must have the 
same form. We find, however, that the material depen- 
dencies of macroscopically quantized fields in linear me- 
dia are not consistent with the classical electromagnetic 
boundary conditions. Subsequently, the boundary condi- 
tions are used to phenomenologically construct quantized 
fields in matter that are expressed in terms of material- 



independent creation and annihilation operators so that 
quantized fields in different media are related by the over- 
all normalization that coincides with the classical bound- 
ary conditions. The generalized position and momentum 
operators of the conforming theory do not obey standard 
commutation relations. We re-derive D'Alembert's prin- 
ciple and Lagrange's equations based on the relativity 
of dynamics in a linear medium [1^, .15.] and show that 
the reduced speed of light changes the definition of the 
canonically conjugate momentum. Then we incorporate 
the modified conjugate momentum into the quantization 
procedure [U, H, H, 0, H, @| and derive quantized fields that 
are the same as the phenomenological fields and likewise 
satisfy the correspondence principle. 



II. QUANTIZATION PROCEDURE 

A typical Ginzburg macroscopic quantization proce- 
dure (6| is based on an expansion of the vector potential 
in terms of modes 



i\ 



(2.1) 



in a dielectric subject to the usual continuum assump- 
tion of an arbitrarily large, linear, isotropic, homoge- 
neous, transparent medium with refractive index n. All 
fields, variables, and operators are macroscopic in the 
sense of representing quantities as continuum averages 
over a volume larger than a cubic wavelength. The 
are orthonormal functions that satisfy periodic boundary 
conditions on the quantization volume. Substituting Eq. 
p.l[) into the electromagnetic Lagrangian 




(2.2) 



and the wave equation, applying separation of variables 
and the identity 

/ (VxU;A)-(Vxu,,A')'^^r= / U,A-[VX(VXU,,A')]'^^I-, 

Jv Jv 

(2.3) 
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one obtains the Lagrangian 

r ^ ( 2-2 2 2 2 \ 



(2.4) 



i\ 



using orthonormality of the mode functions. The conju- 
gate momenta [1, [3, [13 



dL , 



(2.5) 



are then used to construct the effective Hamiltonian 
1 (P'lx 



IX 



,2, .22 



(2.6) 



The effective Hamiltonian takes the form 

^ = ^ E (^'a + ^?Ql) (2.7) 



by apphcation of the canonical transformation 

Pix = Pix/n (2.8a) 



Qix = nqix- 



(2.8b) 



Quantization is accomplished by treating Pix and Qix 
as operators satisfying the commutation relations 



[qix,Pi'x][Qix,Pi'x] = ifi5ii>xx' 



(2.9) 



One then introduces the creation operators and the 
annihilation operators Qx by 

Pix^i\[^{C,l-Qx) (2.10a) 



QiA = ^^(C/A + OA), (2.10b) 

where the operators obey boson commutation relations 

[Oa,0'A'] - [ClxXvx'] = ; [QxXlx'] = SwSxx'. (2.11) 
The effective, or macroscopic, Hamiltonian is therefore 

^=^E^^'(^/aOa+C/a4). (2.12) 
;a 

It is then straightforward to derive the traveling-wave 
representation of the medium-assisted vector potential 
operator [l|, 0, i, 1 i, i, 0, S i, M 



IX 
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The effective Hamiltonian p.l2p and the macroscopic 
vector potential operator (|2.13p are the principle prod- 
ucts of the quantization procedures that appear in the 
macroscopic quantum electrodynamics literature [U, 0, S 

aaiaaasfl]. 

The normalization of the expansion of the vector po- 
tential in Eq. (j2.ip is arbitrary. It is sometimes custom- 
ary to use a normalization of 1/n or 1/ ^/e in Eq. (j2.ip . 
This is equivalent to the canonical transformation (|2.8[) 
and eliminates a step in the derivation. No normaliza- 
tion was used here in order to facilitate comparison with 
a similar quantization procedure that is presented later. 



III. BOUNDARY CONDITIONS 

In the limit of large numbers, quantum and classical 
relations must have the same form. The nature of macro- 
scopic fields as averages means that the limit will be sat- 
isfied. Classically, the spatial extent of a field inside a 
material is reduced in width by a factor of n compared 
to the width in vacuum due to the reduced speed at which 
light travels in the medium. The refracted field is then re- 
duced in amplitude by ^Jn compared to the incident field 
from vacuum as a requirement of conservation of electro- 
magnetic energy in the limit that reflection is negligible 
or suppressed by an anti-reflection coating. 

We consider an electromagnetic field normally incident 
on a minimally reflective interface between the vacuum 
and a linear medium. Forward propagating plane-wave 
solutions of the wave equation can be represented by the 
vector potential 



Af ~ Af cos {—ujt + kz 



(3.1) 



Here, is a unit vector transverse to the direction of 
propagation and k = nuj / c. The amplitudes of the vector 
potential for the incident and refracted fields are respec- 
tively denoted as Ai and At- 

The Poynting-Umov vector cE x B is the continuous 
energy flux vector that is associated with energy conser- 
vation. In terms of vector potential amplitudes, 



(3.2) 



is obtained from continuity of the Poynting-Umov vector 
in the plane- wave cw limit. We write 



At 



Ay 



(3.3) 



for the relation between the refracted field At and the 
field Ai — Ay that is incident from the vacuum. Be- 
cause fields can be transferred from one material to 
another through index-matching or other anti-reflection 
techniques. 



mAi = y/^A2 (3.4) 

is the relation between vector potential amplitudes in dif- 
(2.13) ferent materials and is equivalent to the Fresnel relation 



3 



in the limit that the change in refractive index is suffi- 
ciently small that reflections are negligible. 

We write the vector potential operator ()2.13p as 




and define macroscopic annihilation and creation opera- 
tors 



ai\ 







l\ 



(3.6a) 



(3.6b) 



By substitution from Eqs. p.6p . the macroscopic vector 
potential operator l|3.5p becomes 



ai^e^'- + ale~^^''^)e^,,. (3.7) 



In this form, the material dependence of the classical 
boundary condition p.4p is satisfied by the normaliza- 
tion prefactor of the field in which case the annihilation 
and creation operators obey material-independent boson 
commutation relations 

[ai\,ai>\'] = [aj^,al^,] = ; [aix,aly] = SwSxx', 

(3.8) 

while the polariton operators obey material-dependent 
commutation relations 

Kix, Ci'y] - [CL Ck'] = ; iOxXly] = nSu'Sxx' . 

(3.9) 

This result contradicts the condition (|2.1ip used in the 
macroscopic Ginzburg quantization procedure in Section 
II. Using the operators (|3.6p . one can show that the com- 
mutator (12.911 should be 



[Qi\,Pi'X'] = [qi\,pi' 



inhS. 



WW 



(3.10) 



based on the polariton commutation relations p.9p . The 
scaling of the commutator in Eq. (|3.10p proves that the 
generalized momentum and position variables are not 
canonically conjugate in the usual sense. 



IV. LAGRANGE'S EQUATIONS IN FILLED 
SPACETIME 



characteristics of relativity and Lagrangian dynamics for 
the case in which the effective speed of light is c/n. 

Consider an inertial reference frame S(t, x,y, z) with 
orthogonal axes x, ?/, and z. Position vectors in S are 
denoted by x = (x, y, z). If a light pulse is emitted from 
the origin at time t = 0, then 



(4.1) 



+ - (-t) =0 



describes wavefronts in the S system. Writing time as a 
spatial coordinate ct/n, the four-vector ["l5i | 



X = [ct/n, x) = {ct/n, x, y, z) 



(4.2) 



represents the position of a point in a dielectric-filled 
four-dimensional spacetime. 

Now consider two inertial reference frames, S(t, x, y, z) 
and S'{t',x',y',z'), in a standard configuration [m. uM. 
[20j in which S' translates at a constant velocity v in the 
direction of the positive x axis and the origins of the two 
systems coincide at time t = t' = 0. If a light pulse is 
emitted from the common origin at time t = 0, then 

{xr+iyr+i^'r-{^t'y = (4.3) 

describes wavefronts in the S' system and Eq. (|4.ip holds 
for wavefronts in S. The material Lorentz transformation 

El 



X = 7(x' -I- vt') 



y = y 



(4.4a) 



(4.4b) 



(4.4c) 



(4.4d) 



is derived by the usual methods [l8|, ll9|, |20| , where 

1 

7 = 



(4.5) 



1 - 



The square of the invariant spatial interval As is [Ijl 

(As)2 = (Ax)^ + {Ayf + (Az)2 - {c/nf{Atf. (4.6) 

Multiplying the preceding equation by —1 and taking the 
square root of the result yields another invariant quantity 



In order to derive the generalized momentum variables 
that are canonically conjugate to the generalized coordi- 
nates in a linear medium, we treat space as being en- 
tirely filled with an isotropic homogeneous continuous 
medium that responds linearly to electromagnetic radi- 
ation. Starting from first principles, we then derive the 



cAt = -v/(Ai)2 - (n/c)2((Aa;)2 + [Ayf + {Azf) 

(4.7) 

from which we obtain the interval of proper time 



dT 



dt 

771 



(4.8) 
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Taking the derivative of the position four-vector (|4.2p 
with respect to the proper time, we obtain the four- 
velocity 



dX dK dt ( dx dv dz 

U = — = — ^c\l, , — , 

dr dt dr \ ' d{ct/n)' d{ct/n)' d{ct/n) 



and the four-momentum 

dx dy dz 



moU = mc 1, 



(4.9) 
(4.10) 



d{ct/n)' d{ct/n)' d{ct/n)^ 

where m = "fmg is the relativistic mass. The four-force 

^ _dF _d¥ dt 
dr dt dr 

is derived in a similar manner. 

In the nonrelativistic limit, 7 = 1, r = t/n, we find the 
three- velocity 



u = nx = c 



dx dy dz 



d(ct/n) ' d(ct/n) ' d(ct/n) 
and the three-momentum 

dx dy dz 



p = riTOX = mc 



d{ct/n)' d{ct/n)' d{ct/n) 



(4.12) 



(4.13) 



in a region of reduced light velocity. 

For a system of particles, the transformation of the 
position vector of the i^^ particle to J independent gen- 
eralized coordinates is 



X., = Xi(r; 91,92, • ■ ■,qj), 



(4.14) 



where r = t/n. Applying the chain rule, we obtain the 
virtual displacement 



j=i 



(4.15) 



and the velocity 



dxj 9xj dqj dxi 

j=i 

of the I*'' particle in the new coordinate system. The 
relation 



c?Uj 9xj 



d{dqj/dT) dqj 



(4.17) 



comes from the derivative of Eq. (|4.16p . Substitution of 
Eq. ((4T7| into the identity 



dx, 



u — ^'"'^ + u - ^ ( '^^^ \ (4 18) 

dr \ ' dqj J dr dqj ' dr \ dqj ) 

results in 

dpi dxi d f d ^ 2\ ^ f ^ 2 



in II I — I — Til 1 1 ■ 

dr dqj dr \d{dqj/dT)2 V dq^ \2 ' 

(4.19) 



by application of the calculus. 

For a system of particles in equilibrium, the virtual 
work of the applied forces vanishes and the virtual 
work on each particle vanishes leading to the principle of 
virtual work 



fi • 5xj 



and D'Alembert's principle 

\^ (f _ d-Pi 
^ [ ' dr 







Sx, ^ 0. 



Defining the kinetic energy of the i particle 
we can write D'Alembert's principle (|4.2ip as 



(4.20) 



(4.21) 



(4.22) 



E 



dT 



dT 
dqj 



^dr \d{dqj /dr) 

using Eqs. (|4.15p and (|4.19p . where 

dx, 
dqj ■ 



5qj = (4.23) 



Q. = Ef' 



(4.24) 



If the generalized forces Qj come from a generalized 
scalar potential function V 16], then we can write the 
Lagrange equations of motion 



d f dL 



dr \d{dqj / dr) 



dL_ 

dqj 



0, 



(4.25) 



where L = T — V is the Lagrangian in a linear medium. 

For any coordinate qj that measures a linear displace- 
ment of a particle in a given direction, the partial deriva- 
tive of the Lagrangian with respect to that coordinate 
vanishes. In that case, substitution of 



dL 

dqj 



= 



into Eq. (|4.25p produces 

d dL 



dr d{dqj/dT) 



= 0. 



(4.26) 



(4.27) 



If the potential energy is velocity independent then we 
can define the canonical momentum 



Pi 



dL 



dL 



d{dqj/dT) c d{dqj / d{ct / n)) 



(4.28) 



based on the proper temporal integration of Eq. (|4.27p . 

The consequence that an effective speed of light has 
for the canonical momentum is easily illustrated in the 
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case of a nonrelativistic free particle. The kinetic energy 
of the particle, initially at rest in the local frame, is 



1 

T = —mu ■ u. 

2 



(4.29) 



Then the Lagrangian for the free particle can be written 
as 



m / dxj 



(4.30) 



in rectangular coordinates. Applying the definition of the 
canonical momentum (j4.28p . one obtains 



dL 



Pj = 



d{dxj/dT) 



dxj dxn 
m— — — mn—- 
dr dt 



—J- = mn—^ (4.31) 



in agreement with the linear momentum of a particle in 
a linear medium, Eq. (|4.13[) . Although macroscopic par- 
ticles cannot travel unimpeded through a material in the 
continuum limit, light does at reduced speed. In the next 
section, the modified Lagrangian dynamics are applied to 
the derivation of macroscopic quantum electrodynamic 
principles. 



V. MACROSCOPIC QUANTIZATION 

Lagrangian dynamics is the basis for the procedure 
that was used in Sec. II to quantize the electromagnetic 
field in a linear medium. The macroscopically quan- 
tized field, Eq. (I2.13p . was found to be inconsistent with 
the classical electromagnetic boundary conditions. In 
the preceding section, we derived a modification of La- 
grangian dynamics in the context of a uniform c/n speed 
of light. Here we show that the quantized field that is 
derived using the modified dynamical theory is consistent 
with the electromagnetic boundary conditions. 

The macroscopic quantization procedure is based on 
an expansion of the vector potential in terms of modes 



A = c^qix{t)uix{r) 



(5.1) 



IX 



Substituting Eq. (|5.ip into the Lagrangian 




(5.2) 



and the wave equation, one obtains the Lagrangian 

^ = ^ E - n'ufql) , (5.3) 

as in Section II. The point of departure for this deriva- 
tion is the use of Eq. rather that Eq. (^31) . for 
the conjugate momenta. Instead of pix = n^qix, which 



leads to violation of quantum-classical correspondence, 
we obtain 



dL 



Pix 



nqix 



(5.4) 



d{dqix/dT) 

for the conjugate momenta. The effective Hamiltonian 
^ = E^^'A^ - ^ = J E (P?A + n^-Ulx) (5.5) 



is quantized in the usual way by taking Pix — \/npix and 
Qix — lixl \/n to be operators satisfying the material- 
independent commutation relations 



XX' 



(5.6) 



Defining the usual material-independent annihilation and 
creation operators, 



1 



a/A 



.{loiQix + iPix) 



4x- 



1 



yields the effective Hamiltonian 



oJiQix - iPi 



IX) 



H 



IX 



l^aix + aixal^ 



(5.7a) 



(5.7b) 



(5.8) 



In the plane- wave cw limit, the classical energy density 
is 

n^-(^[--) (V X A)^ j = -^A^. (5.9) 

Using the classical boundary conditions At — Ay/^/n, 
the energy density inside the material is a factor of n 
greater than the energy density in the vacuum. The 
effective Hamiltonian (|5.8p exhibits the conforming en- 
hancement of energy density in a medium. In contrast, 
the effective Hamiltonian (|2.12p in the material, derived 
by the original Ginzburg quantization procedure with the 
boson commutation relations, fails the test of quantum- 
classical correspondence. 

The traveling wave representation of the macroscopic 
vector potential 



A = cE 



(A 



2nLJiV 



aixe 



+ at,e-'^'-)eu,, (5.10) 



can be constructed from the relations (|5.7p between the 
canonical variables and the creation and annihilation op- 
erators using the canonical transformation. The field 
(I5.10p . like the phenomenological field p.7p . satisfies 
quantum-classical correspondence with the electromag- 
netic boundary condition At = Ay/y^. 
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VI. SUMMARY 

The quantized vector potential and the effective Hamil- 
tonian for an electromagnetic field in a dielectric were 
derived from the Lagrangian by the conventional quan- 
tization procedure. This representation of the macro- 
scopically quantized field was shown to violate quantum- 
classical correspondence. The principles of Lagrangian 
dynamics, including D'Alembert's principle and the La- 



grange equations of motion were re-derived using the 
thesis of continuum electrodynamics of a structureless 
medium that interacts with light through a constant. 
The resulting changes to the principles dynamics in the 
continuum, applied in the Ginzburg macroscopic quanti- 
zation procedure, were found to be sufficient to repair the 
violation of quantum-classical correspondence for macro- 
scopically quantized fields. 
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